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Abstract. Let V be a set of isolated points in R''. Define a Unear functional £ 
on the space of real polynomials restricted on V, Cf = '}2xev /(^)p(^)> where 
p is a nonzero function on V. Polynomial subspaces that contain discrete 
orthogonal polynomials with respect to the bilinear form (/, g) = C{fg) are 
identified. One result shows that the discrete orthogonal polynomials still 
satisfy a three-term relation and Favard's theorem holds in this general setting. 



Discrete orthogonal polynomials appear naturally in combinatorics, genetics, sta- 
tistics and various areas in applied mathematics (see, for example, In one 
variable they have been studied extensively. Let y be a set of isolated points on 
the real line, its cardinality \V\ is either finite or countable. Let w be a real positive 
function on V. With respect to the bilinear form (/, 5) = X^^ev 
there is a sequence of orthogonal polynomials {p„ : < n < \V\} on V with 
{PmPm) = for n ^ m (for example, using Gram-Schmidt process). These are the 
discrete orthogonal polynomials. Their structure is similar to that of the usual con- 
tinuous orthogonal polynomials. For example, every sequence of discrete orthogonal 
polynomials satisfies a three-term relation, 

(1.1) Xpn^ anPn+l+KPn + CnPn-l, 0<n<|y|-l, 

where a„, 5„ and c„ are real numbers. Furthermore, according to Favard's theorem, 
the three-term relation essentially characterizes the orthogonality of polynomials. 

Discrete orthogonal polynomials of several variables are far less studied. Their 
orthogonal structure is much more complicated than that of one variable. Even 
some basic problems have not been addressed. Let us first fix some notation. We use 
the standard niultiindex notation: for a; G M'^ and a E Wq, write 
This is a monomial of (total) degree \a\ := ai -I- . . . -I- ad- Let li^ = R[xi, . . . , Xd\ 
and let IV^ be the subspace of polynomials of degree at most n. Denote by the 
space of homogeneous polynomials of degree n. It is well known that 



Let £ be a linear functional defined on 11 , such that a basis of orthogonal poly- 
nomials {Pa : |q;| = n, a G No,n > 0}, where Pa G Ilf^, exists with respect to 
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the bilinear form {f,g) ~ C{fg), in the sense that {Pa,Pf3) = if \a\ ^ \(3\. Let 
V,^ = spanjPa : \a\ — n} be the space of orthogonal polynomials of total degree n. 
Then dim V,'J ~ dim'P,'^. In this case, if we adopt the point of view that the orthog- 
onality holds in terms of the subspaces V^, not in terms of particular bases of V^, 
then we can have an analog of a three-term relation. Let P„ = {Pq, : \a\ — n}; we 
also use P„ to denote a column vector, in which the elements are ordered according 
to a fixed monomial order. Then the following three-term relation holds, 

+ B„,,P„ + Cna^n-l, 1 < * < > 0, 

where An^i, Bn,i and Cn,i are matrices of appropriate dimensions, and P_i :~ 0. 
Furthermore, there is an analogue of Favard's theorem ( 1 1 ) ■ For the general theory 
of orthogonal polynomials of several variables, we refer to Chapter 3 of 0. 

Let y be a set of isolated points in M"*. Again we denote by |y| the cardinality 
of V, which can be finite or countable. The orthogonal polynomials on V depends 
on the structure of the polynomial ideal I{V) that has V as its variety, 

liV) ^{pe M[xi, ...,Xd]: p{x) = 0, G V}. 

The discrete orthogonal polynomials on V can only consist of polynomials that 
do not belong to I{V). Let C be defined by £f — J2xgv where W is 

a real function on V, W{x) ^ for aU x £ V, and 2a;ey < oo for 

all a S Nq in the case where is a countable set. Only if I{V) = {0}, can a 
complete basis {Pa ■ a € Ng} of IT^ with respect to £ exist and the discussion in 
the previous paragraph apply. If the ideal I{V) is non-trivial, for example, when 
\V\ is finite, then we need to understand the subspace M.[V] = Il'^/I(V) in order 
to define orthogonal polynomials. In the case M.[V] = 11^, for example, little extra 
work is needed; the three-term relation in the form (|1.2(l holds for < rt < A^. This 
is the case of straightforward extension of discrete orthogonal polynomials in one 
variable. However, even in the case where is a product of two point sets X and 
Y in one variable, V = X x Y with \V\ = N, M.[V] consists of only a subspace of 
n^. In general, the space R[V] can be rather complicated and care is needed for 
the definition of orthogonal polynomials on a discrete set V. 

The purpose of the present study is to define discrete orthogonal polynomials in 
this general setting. The polynomial subspaces for which the discrete orthogonal 
polynomials exist are identified and the three-term relation and Favard's theorem 
are established. In the following section we discuss the structure of polynomial 
subspace on V. Discrete orthogonal polynomials on V are studied in Section 3. 
Various examples will be given in the paper, some are given in terms of the classical 
discrete orthogonal polynomials, such as Hahn polynomials, in Section 4. It is our 
hope that this study can help to clarify some of the basic questions in the theory 
of discrete orthogonal polynomials. 

The author would like to thank Paul Terwilliger for the discussion on the Leonard 
pairs and tridiagonal pairs (see |7| |Sj and the reference there), during the 7th 
International Symposium on Orthogonal Polynomials and Special Functions, held 
in Copenhagen in August, 2003, which suggests this study. 

2. Polynomial spaces on V 

First we need to understand the structure of the quotient ideal 11''//, where 
/ :— I{V) and F is a set of isolated points in K*^, finite or countable. Most of the 
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results below also hold if V has finitely many accumulation points. We review some 
results about ideals and varieties, our basic reference is 

For f,g e n'', we say that / is congruent to g modulo /, written as f = g mod /, 
if and only if f — g € I. If |y| is finite, then it is known that the codimension of 
/ is equal to \V\; that is, dimW^ /I{V) = \V\. Let R[V] denote the collection of 
polynomial functions cf) : V t-^ R. This is a commutative ring and it is isomorphic to 
the quotient ring Il'^/I{V). We can identify M.[V] with 11''// as there is an one-to-one 
correspondence between (p S M.[V] and [(p] ~ {g ^ H"^ : g = (p mod /}. It is possible 
to say more about this space. For a fixed monomial order, we denote by LT(/) the 
leading monomial term for any polynomial / G H"^; that is, if / = ^CaX°', then 
LT(/) — Cj3X^ , where is the leading monomial among all monomials appearing 
in / and cp ^ 0. For a polynomial ideal / other than {!}, we denote by LT(/) the 
leading terms of /, that is, 

LT(/) = {cx"|there exists / e / with LT(/) = cx"}. 

We further denote by (LT(/)} the ideal generated by the leading terms of LT(/) for 
all / G / \ {0}. According to the Hilbert basis theorem, every polynomial ideal has 
a finite basis. A set {gi, . . . , } is called a Grobner basis of / if 

(LT(5i),...,LT(gO) = (LT(/)). 

It is known that every polynomial ideal has a Grobner basis. Such a basis enjoys 
many interesting properties that have important applications. For example, it is 
used to prove the following result Chapt. 5]). 

Proposition 2.1. Fix a monomial ordering on IV^ and let I C 11'' he an ideal. 
Then there is an isomorphism between 11''// and the space 

Si :=span{x"|x" ^ (LT(/))}. 

More precisely, every f G is congruent modulo I to a unique polynomial r G 5/ . 

In fact, the polynomial r is the remainder of / on division by /. For / = 
and a fixed monomial order, the division algorithm states that for 
every / G 11'', there exist pi and r in 11'' such that f = J2 Pifi + '^'j where r G 5/ and 
no term of r is divisible by any of LT(/i), . . . , LT(/Af). The remainder polynomial 
r is unique if the basis /i, . . . , is a Grobner basis. 

Proposition 2.2. Let V be a set of isolated points in R"^ and I = I{V). Let 
A :~ A(T^) be the index set A — {a : x"^ ^ LT(/)}. Then every polynomial P G M.[V] 
can be written uniquely as 

P{x) = ^Ccx" mod /(F), CaGR, 

aeA 

and the set A satisfies the following property 

(2.1) a G A implies a — (3 £ A, whenever a — G Ng and f3 G Nq. 

Proof. For / = I{V), we can take Si in Proposition 12.11 as R[V], modulus / if 
needed. The definition shows that we can write Si = spanja;" : a G A}. Hence, 
every polynomial P in M[y] has the stated representation. Since the ideal (LT(/)) 
is a monomial ideal, x" G {LT(/)) implies x"~^^ G (LT(/)) for any (3 G Nq. Conse- 
quently, it follows that the set A satisfies the property (j2.1|l . □ 
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In the following, we shall drop modulus / and use M.[V] to denote the space 

(2.2) R[V] = span{x" : a G A(F)}. 

This abuse of notation should not cause problems. 

We should point out that the set A is not unique, since all equations actually 
hold under congruence modulo /. In fact, Grobner bases are not unique, since the 
choice of monomial orders matters. There are in fact many different representations 
of elements in M.[V]. What is of interest is the property (|2.1I) satisfied by A. 

Example 2.1 Consider the set V = {(0,0), (0, 1), (1,2), (2,3)}. It is easy to see 
that I{V) = (51,52,53), where 

gi{x,y) ^ x{x -l){x -2), g2{x,y) ^ x{x + 1 - y), g3{x,y) ^ y{x + 1 - y). 

If we use the graded reverse lexicographical order (n — fc— 1,/c+l) >- {n — k,k), 
then (LT(J)) = {y'^,xy,x^) and R[V] — spa,n{l,x,y,x^}. If we use the graded 
lexicographical order {n — k,k) >- (n — fc — l,fc + 1), then (LT(/)) = {x'^,xy,y^) 
since we also have I{V) = (52,53,54) where gi{x,y) — y^ — 6y^ + 5y + 6x, and 
M[y] = span{l,a;,y,?/2}. □ 

For d — 2, the property H2.1|l of A shows that the set A must be of a stair shape 
as the lattice points in the unshaded area depicted in Fig. 1. 

8' 1 1 1 1 1 1 1 1 =1 1 1 1 1 1 
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2 4 6 8 1 2 3 4 5 

Figure 1 

That is, in the case of d = 2, for each set A there is a sequence of positive integers 
rii, which satisfies n„j < rim-i < • ■ • < (some of the Ui can be positive infinity), 
such that 

(2.3) A = {(fc,0 : < / < m,0 < A: < n,}. 

Example 2.2. Let A be the lattice set in the two figures in Figure 1. For the 
left one, m = 6 and A = {(i,j),0 < i < nj,0 < J < 6} with (no,...,n6) = 
(6, 4, 4, 4, 2, 2, 0). For the right one, to = 3 and 

A = {(i, 0) : < i < 2} U {{i, 1) : < i < 2} U {(0, 2)} U {(0, 3)} 

with (no. Til, 712, nj) = (2, 2, 0, 0). □ 

Proposition 2.3. There exists a point set V for which M.[V] in H2.2|l is given by 
the index set A in H2.3() . 
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Proof. Let xq, xi, . . . , Xng and j/o, ■ • ■ , 2/m be isolated real numbers. We define 
the point set V as follows: 

V ^ {{xo,yo), {xi,yo), (x„o,?/o), 
{xo,yi), {xi,yi), {xni,yi) 

{xQ,ym), {Xi,ym), 

If no is finite, then \V\ ~ no + ni + . . . + n„i + m+l. We let Um+i — and adopt 
the convention that Y[i=o '^i = 1 • Define polynomials 

5fe(a;) = - a;i) ]^(?/ - yj), < k <m + l, 

i=0 j=0 

where if Uk = oo, then we take gk{x) = 1. Then it is easy to see that (/) = 
(50,51, ■ • ■ ,5m+i), so that 

(LT(/)) = . . . , 

from which it follows that M.\V] = {x^y'- : (k, I) G A} with A given in p. 3(1 . □ 

Example 2.3. Consider the "triangle" point set 

y^{{xQ,yo), {xo,yi), (a;o,?/m); 

(xi,yi) {xi,ym) 

where all points are isolated and |y| = (m + l)(m + 2)/2. In this case, ni — m — I 
for < / < TO (making a transpose of the array as in matrix transpose), so that 
A = {(fc, : < fc + / < m} and Wy] = H^^. □ 

Example 2.4. Consider the "product" point set 

V = X Y.Y = {(xj,yj) : < i < n,0 < j < m}, 

for which \V\ = (n + 1)(to + 1). In this case = ni = . . . = rim = n- and 
R[V] = {x'^y^ ■.0<k<n,0<l<m} = nlx H^. □ 

The above discussion can be extended to M[a;i, . . . ,Xd]- We stick to the case 
d = 2 to keep the notation simple. 

A couple of remarks are in order. First, if to = cxd in Example 2.2, then K[V"] = 
R[x,?/]. This is also the case for Example 2.3 when both n and m are infinity. If, 
however, n is infinite and to is finite (or other way round), then A = {(^,0 • ^ 
Z < m, fc > 0} is an infinite set but not all Nq and R[V] = 11'' x 11™. Furthermore, 
the space of polynomials of one variable appears as a special case. 

Example 2.5. If T/^ is a point set on the x coordinate line, V — {{xq, 0), . . . , (a;„, 0)} 
{n can be infinity), then R[V] = 11}^. □ 

Second, we should point out that the point set V in the proposition, and the 
above examples, are the simplest examples for which M.[V] can be determined. For 
a generic point set V, the problem of determining is highly non-trivial. One 
possible algorithm, at least for V is finite and moderate in size, is to check the rank 
of the matrices whose rows are the vectors Xa ■— {x" : x G V}. 
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The algorithm goes as follows: Fix a graded monomial order, starting with 
a S Nm '■— {/? : \(3\ < m} for a small m, so that the resulting matrix {Xa)aemrn 
has full rank. Then add new rows Xa according to the order in |q!| = m + 1 to the 
matrix. For each new row added, check the rank of the new matrix; if it has full 
rank, add the next Xa and proceed; if it does not have full rank, remove this row 
and add the next X^ and continue. When the matrix becomes a square nonsingular 
matrix, the corresponding set of x" will be a basis for R[F]. 

3. Discrete orthogonal polynomials 

Let F be a set of isolated points in K'^. Let be a real function on V and 
W{x) ^ for any x ^ V. Assume that 

^ |x"||M^(a;)| < oo for aU a e NfJ 

xev 

in the case where V is an infinite set. We define a bilinear form (•, •) on H'' x H'' by 

(/,5)=A/5), where ^ /(a;)W^(.T). 

xev 

If {f,g) = 0, we say that / and g are orthogonal to each other with respect to 
W on the discrete set V. The notation C{fg) is more convenient for the matrix 
operations below. 

Fix a graded monomial order. Let M.[V] and A — A{V) be defined as in the 
previous section (see H2.2ll ^. Let n — max{|a| : a S A(y)} and Ak{V) — {a € 
A{V) : |a| = fc}, < fc < n. Note that n can be infinity. Define = |Afc(F)|. We 
denote by x'^, and x*^, < fc < n, the sets 

x^ = {.t" : a e AiV)}, x^ = {x" £ x'^ : |a| < fc}, x'^ = {.t" : a £ AkiV)}, 

respectively. We also regard them as column vectors in which the elements are 
ordered according to the fixed graded monomial order. 

The set of orthogonal polynomials on V will be denoted by {Pa, a £ A{V)}, 
where Pa has degree \a\. We introduce the following notion. If {Pa : a £ A{V)} 
is a sequence of polynomials in R[V], then set Pfe := {Pa ■ a £ Ak{V)}. Just as in 
the case of x'^, we also regard as a column vector. 

Definition 3.1. Let V be the set of isolated points and W be a nonzero real function 
on V as above. A sequence of polynomials {Pa £ n|Q,|d : a £ A(V^)} is orthogonal 
with respect to W if 

£(x'P^) = (x',Pfc) =0, k>l, and £(x''Tj) = (x^ P^) = 5^, 

for < k < n, where Sk is an invertible matrix of size r^ x r^. The sequence 
is orthonormal with respect to W on V if (Pfc,Pfc) = Ij.^, the identity matrix, for 
< fc < n. 

The notation £(x'P|") is more convenient than (x'^,Pfc}, since it shows clearly 
that this is a matrix of size rfe x rfe. The orthogonality of Pa is defined as orthogonal 
to lower degree polynomials, as in the continuous case. The polynomials of the same 
degree may not be pairwise orthogonal. By definition, we can write 

Pfc = Gfex'^ + Gfc_i,fcx'=-i + . . . + Gi,fex°, 

where Gk is a x matrix, called the leading coefficient of P/j. Assume a sequence 
of orthogonal polynomials Pa exists on V . Then we can follow the proof in |21 
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Section 3.1] to show that {Po, . . . ,Pfc} is a basis for K[V^] nn^, Hk := C{PkPl) and 
Gk are both invertible. Furthermore, the following theorem still holds. 

Theorem 3.2. A sequence of orthogonal polynomials {Pa '■ a G A{V)} with respect 
to W on V exists if and only if the matrices Mk '■— (x^,x^) are nonsingular for 
< fc < n. 

The proof is based on linear algebra and follows exactly as in the continuous 
case; see |2 Theorem 3.1.6]. 

If W is positive on V , W{x) > for all x € V, then (/, f)—0 implies / = on 
V, so that the bilinear form (•,•)= becomes an inner product on M.[V]. For such 
a W, orthogonal polynomials on V exist. Furthermore, in this case, we can have 
orthonormal bases. 

Theorem 3.3. If W is a positive function on V , then a sequence of orthonormal 
polynomials {Pa '■ oi G A(y)} with respect to W on V exists. 

Proof. In this case, the matrix Mj. is positive definite since for any nonzero column 
vector c, MkC = (c-^x'^, c^x*^} > 0. In particular. Ma := (x^, x"^) is a symmetric 
and positive definite matrix. It follows that it can be factored as Ma = SDS"^ where 
5 is a nonsingular lower triangular matrix and D = diag{(ii, ^2, . . . , c?|y|} with all 
d, > 0. Let D^i/^ ^ diag{^^~^/^ . . . , d~yl^} and R = D-^/^S~^. Then RM^R^ = 
{Rx.^,Rx^) = Since S is lower triangular and x'^ = {x*', x^, . . . , x"} as a 

column vector, we can write the components of i?x^ as Pq, . . . , Pn where Pfc consists 
of polynomials of degree k. These are the orthonormal polynomials. □ 

The proof of the theorem provides an algorithm that can be used to construct 
discrete polynomials in several variables. For \V\ is finite and moderate in size, it 
is rather effective; an example is given in the following section (Example 4.3). 

Let Vk{W) denote the space of orthogonal polynomials of degree k; that is, 
Vk{W) = spanPfc. Evidently dimV/c(VF) = r^. Comparing to the orthogonal 
polynomials in the continuous case, the numbers depend on V and there is no 
closed formula for them. Furthermore, r/j < dimP^ and the equality often does 
not hold. In fact, may no longer be an increasing sequence. 

Example 3.1. Let V = Xy.Y \>c the product point set in Example 2.3, where X — 
{xo, . . . ,a;„} and Y = {j/o, ■ • ■ , 2/m}- Here n = n + m. Let W(x,y) = Wi(a;)u;2(y), 
where W\ is positive on X and W2 is positive on Y . Then the orthogonal polynomials 
Pfc, < fc < n exist and can be constructed as follows: Let {poiPi, ■ • ■ jPn} and 
{go, 91, • ■ • , qn\ denote the discrete orthogonal polynomials with respect to w\ on 
X and W2 on y, respectively. Then the orthogonal polynomials on V are given by 
Pk,i{x,y) = Pk{x)qi{y). We assume m is finite and n > m. Then in our vector 
notation, 

^k{x,y) = {po{x)qk{y),pi{x)qk-i{y), . . . ,Pk{x)qo{y)} for < fc < m, 
so that rfc = fc + 1 for < fc < m, 

Pk{x,y) = {Pk-m{x)q„i{y), . . . ,Pkix)qoiy)} for m + 1 < fc < n, 
so that rk — m + 1 for m + 1 < fc < n, and 

Pk{x,y) = {Pk-rn{x)q„i{y), . ■ . ,Pn{x)qk-n{y)} ioT n + 1 < k < n + m, 
so that rk — n-\-m — k + 1 for n + l<k<n + m— 1. □ 
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Next we consider the three-term relations satisfied by the orthogonal polyno- 
mials. If Pa is an orthogonal polynomial, then Pa G M.[V] so that it is a linear 
combination of x" for a € A(V). Clearly, multiplying by a coordinate Xi gives a 
polynomial XiPa of degree \a\ + 1. However, unlike the continuous case, Xix" may 
not belong to M.[V] for some a £ A(y). Nevertheless, it is congruent modulus I{V) 
to a unique polynomial in K.[V]. Recall that n = max{|Q!| : a G A(y)}. 

Proposition 3.4. Let I := I{V). For < fc < n — 1, there exist matrices Ak^i : 
rk X rfc+i, Bk^i ■■ rk xrk, and Ck,i : rk x rk-i, such that for 1 < i < d, 

(3.1) XiPfc(x) = Ak,iFk+i{x) + Bk,^Fk{x) + C\,Pk-i{x) mod /, 

where < fc < n — 1 and we define P_i = 0, Aks = and C-i^i = 0; moreover, 

(3.2) Ak,^Hk+i = /:(.T,;PfePj+i) = FfeCj+i^„ S,,,Ffe = Cix.lPkFl). 

Proof. If all components of .x^Pfe arc in K[y], this is proved as in the usual case, 
by writing XiFk in terms of orthogonal polynomials Po,Pi, . . . ,Pfc+i and then use 
orthogonality. If Pa £ M.[V] but XiPa ^ M.[V] for some |a| = k, then there exist a 
Q in I{V) and an Ra £ M.[V] such that XiPa{x) = Q{x) +Ra{x), and the degree of 
Ra is at most k + 1. We can write Ra = ak+iPk+i + afePfe -I- ak-i¥k-i + ■ ■ ■, where 
Bj are some row vectors of appropriate size. By the orthogonality, we get 

ak+lC{Pk+lPl+l) = C{Ra¥l+^) = C{PaPl+,), 

since Q vanishes on V and (Q, P) = for any P. Similarly, we get afeiJfe = 
£(PaP^), ak-iHk-i = £(PqP|^_j), and all other aj are equal to zero. In vector 
and matrix notation, this is the three-term relation. The presence of Q means that 
the equality holds under modulus I{V) in general. □ 

Corollary 3.5. Let I := /(V). // {Pfc} o,re orthonormal polynomials, then for 
l<i<d, 



(3.3) Xi¥k{x) = 

where < fc < n 
symmetric. 



^fc,^Pfe+i(a;) + Bk,^Pk{x) + Al_,^^Fk-i{x) mod /, 
- 1, P_i = 0, Ak,i = and C-i^i = 0; moreover, Bn^i are 



In the case of continuous orthogonal polynomials, the matrix A„ ^ has more 
columns than rows and it has full rank. This is no longer true in the discrete case. 
Since rk may no longer be an increasing sequence, the matrix An^i can have more 
rows than columns; moreover, it may not have full rank. 

Example 3.1*. We continue the example V ^ X xY in Example 3.1. Assume 
that {po,Pi, ■ ■ • ,Pn} satisfies the three-term relation 

xpkix) = anPk+iix) + bkPk{x) + CkPk-i{x) 
and {go) 9i) • • • j 9m} satisfies the three-term relation 

yikiy) = a^qk+i(y) + b'kqkiv) + c'^qk-iiy), 
respectively. With x — xi and y = X2, the matrices Ak^i and Ak^2 take the form 



A 



fc,i 



ao 



O 



01 



dk 



O 0] 



and Ak,2 



0<k <m, 



o 
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of size (fc + 1) X (fc + 2); 



a-k- 



o 



O' 



and Ak,2 = 



0' 

o 



o 



< fc < 



of size (m + 1) x (m - 

^k—rn 



A. 



k.l 



1); and 

01 



o 





Ofe 





and Ak.2 





o 



o 



n < k < n + m, 



of size {n + m — k + 1) x (n + m — k). Note that Ak.2 in the case m < k < n does 
not have full rank. For n < k < n + m, Ak^i has more rows than columns. □ 

One consequence of the three-term relation is the Christoffel-Darboux formula, 
h+i{x)Al ^Fk{y) - Pk{x)Ak,^Pk+l{y) 



k 
3=0 



'iix)Fkiy) = 



mod I{V), 



where 1 < i < d and < fc < n. The proof follows as in the continuous case. 

The composite matrix Ak — (^^"1, ■ • ■ , ^fe d)^ plays an important role in Favard's 
theorem of several variables. This matrix is of size dvk x r^^+i. 

Proposition 3.6. For < fc < n — 1, drk > rk+i; the composite matrix Ak of 
Ak^i, . . . , Ak.d and the composite matrix Ck+i of Ck+i,i, ■ ■ ■ , Ck+i.d both have full 
rank, 



(3.4) 



ranky4fe — rankCj_^]^ = fk+i- 



Proof. Recall that Gk denotes the leading coefficient of and it is an invertible 
matrix. Let a = (0, . . . , 0, 1, 0, . . . , 0) be the «-th element of the standard Euclidean 
basis. By a € A{V) implies that a - et e A(y). If a e Ak{V) but a + ^ 

Ak+i{V), then Xix"" = 2;"+'^* £ (LT(/)). We then define the matrix Lk,i by 

x,x'' = Lfe^.x'^+i mod (LT(/)), 1 < i < d. 

The matrix Lk^i is of the size x rk+i and it is uniquely determined. Comparing 
the coefficients of x'"'+^ in both sides of the three-term relation, we see that 



(3.5) 



GhLk.i — Ak.iG 



kA'^k+li 



l<i< d. 



The entries of the matrix Lk^i are mostly with at most one 1 in each row. However, 
Lk,i may not have full rank. 

Since a € Ak+i{V) implies that a — Ci € Ak whenever — 1 > 0, it follows that 
drk > rk+i since — \Ak{V)\. Moreover, since the column vector (xix'^, . . . , Xdx'^) 
is equal to Lfex''+^, where Lk is the composite matrix of Lk,i, ■ ■ ■ , Lk,d, and clearly 
{xiX*^ : 1 < i < d} includes x'^'^^ as a subset, it follows that Lk has full rank rk+i- 

The equation H3.5|l imphes that AkGk+i = diagjGfe, . . . , Gk}Lk- Since Gk in- 
vertible implies diag{Gfe, . . . , Gk} invertible, it follows that rank = rk+i- Fur- 
thermore, H3.2|l implies that AkHk+i = diagjiJ^, . . . , Hk}G^^^ and Hk is invertible; 
hence, rankCfe+i — rankA^. □ 
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Since the matrix Ak has fuU rank and drk > rk+i, it has a generahzed inverse, 
, which is of the size rk+i x drk and can be assumed to be of the form = 
[DJ.^, . . . , DJ^ where - are of the size rk+i x r^. Then 

D^Ak = Dl^Ak,i + . . . + Dl^Ak,d = /r.+i- 

We note that the generahzed inverse is in general not unique. Using D^, we get 
from the three-term relation a recursive formula 

d d d 

(3.6) Pfe+i = ^^'^k,^^k - ^k'^^f^ - Ck,^Pk-l, 

i—1 i—1 i—1 

which allows us to compute Pfc+i using P^ and Pfc-i. This formula is useful in the 
proof of the analog of Favard's theorem. 

According to Propositions 3.4 and 3.6, orthogonal polynomials on a set V of iso- 
lated points satisfy a three-term relation whose coefficient satisfies a rank condition. 
We want to establish that the converse is also true, that is, an analog of Favard's 
theorem. To this end, we start with a sequence of polynomials that satisfies the 
three-term relation H3.1|l and the rank condition (|3.4I) . and show that there exist a 
set V of isolated points and a weight function on with respect to which the 
polynomials are orthogonal. 

For this purpose let us start with an ideal / C 11'' and let A := {a G Nq : 
x" ^ (LT(/))}. The proof of Proposition 12.21 shows that A satisfies the property 
(12. Assume that there is a sequence of polynomials Pa £ 5/ := span{x"|x" ^ 
(LT(/))}, where Pa is indexed by a e A such that Pq — 1 and Pa G nj^^|. Set 
n = max{|a| : a G A}, which can be infinity. Let A^ = {a G A : |q:| = k} and let 
Pfc = {Pq, : a G Afc} for < fc < n, and regard Pfc as column vectors according to 
a fixed graded monomial order. 

Theorem 3.7. Let I be an ideal o/II'^ and let A and P^ be as above. Assume that 
Pfc satisfies the three-term relation (13. 1() whose coefficient matrices satisfy the rank 
condition (|3.4|l . 

(i) There is a linear functional C on Si for which Pa are orthogonal polynomials 
with respect to the bilinear form C{fg) — {f,g)- 

(ii) If n is finite then there exist a set V of isolated points and a real function W 
on V, W{x) ^ for all x € V , such that {Pa : a G A} is a sequence of orthogonal 
polynomials with respect to W on V . 

Proof, (i) The proof follows along the line of the proof of Favard's theorem for 
continuous orthogonal polynomials of several variables. We shall be brief whenever 
proofs in the two cases are essentially the same. Using induction, it follows from 
the three-term relation and the rank condition that the leading coefficients Gk of 
Pfc are invertible. The linear functional, £, defined by 

£1 = 1 and C{Fk) ==0, 1 < fc < n, 

is well-defined for Si, since every polynomial in Si takes the form X^aeA ^ax" and A 
satisfies H2.1II . Furthermore, using 1)3. 6|l . one can show by induction that C satisfies 
£(PfcPj) = for fc ^ j, and the matrix Hk = £(PfcP^) is invertible. Consequently, 
Pfc are orthogonal polynomials with respect to the bihnear form C{fg) — {f,g). 
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(ii) We only need to show that the hnear functional C can be represented by a sum 
over a set of isolated points; that is, C can be written as 



for some with Ai 7^ for 1 < i < A^. Assume that C has such an expression. 
If Xi are known then CPa = ^a,Oj Q? 6 A, becomes a system of equations on A^. 
The coefficient matrix of this system is [Pc((xi)], where a £ A and 1 < i < A^. 
In particular, if A'^ = |A| then the matrix is a square matrix. Its determinant is a 
polynomial in variables Xi,X2, . . . ,xjv and defines a hypersurface in R'^^. Hence, 
for almost all choices of the values of Xi,...,xjv, the determinant is non-zero. 
Furthermore, by Cremer's rule, it is possible to choose a set F = {xi, . . . , xat} such 
that Xi ^ Q ioT 1 < i < N . Hence, Pa are orthogonal with respect to W on V , 
where W is defined by VF(xi) = A,, 1 < i < iV. □ 

A couple of remarks are in order. First of all, if n is infinity, we do not know if 
C can be written as a sum of point evaluations, Cf = X^xgv /(■^)' '^'^ ^ countable 
set V . In fact, if / is the trivial ideal (1), then A — Nq; the three-term relation 
takes the same form as that for continuous orthogonal polynomials and the rank 
condition remains the same. Hence, the three-term relation and the rank condition 
are not enough to give further information on the linear functional. 

The same phenomenon will happen to the case that Si is a product polynomial 
space, say, HJ^ x H'', for which the orthogonal polynomials are pj {x)qk (y) , < j < m 
and fc > 0. No matter if qk are orthogonal with respect to a linear functional defined 
by an integral or to a linear functional defined by an infinite sum, the three-term 
relation will take the same form and the rank condition will also remain the same. 

Secondly, the proof of the theorem only shows that W is nonzero at every point 
of V. This is enough if we only deal with orthogonality but not orthonormality. 
The theorem simply states that the rank condition H3.4() and the three-term relation 
(|3.1() are enough to ensure orthogonality. The following corollary is about the case 
of orthonormality, where we do get positive weight. The difference is in the three- 
term relations (|3.3|l vs 

Corollary 3.8. Let I be an ideal of H"^ and let A and be as in the theorem. 
Assume that n is finite. IfPk satisfies the three-term relation ()3.3() whose coefficient 
matrices satisfy the rank condition (13.4(1 . then there is a set V of isolated points and 
a positive function W on V such that {Pa : a G A} is a sequence of orthonormal 
polynomials with respect to W on V . 

Proof. According to the theorem, there is a set V so that Pa are orthogonal with 
respect to the bilinear form defined by the linear functional C of the form 1(3. 7|l . 
We need to show that Hk = £(PfcFj) is an identity matrix for < fc < n. This 
can be established by induction. Since Pq = 1 and £1 = 1, we have Hq = 1. By 
(E3 with Cj+i = Ak, AkHk+i = diag{iJfc, . . . , Hk}Ak. Assume = U^. Then 
diag{iJ/c, . . . , -fffe} is an identity matrix, so that Hk+i is an identity matrix by the 
rank condition. From the fact that £(PfeP|") = /, it follows easily that £ is a 
positive definite linear functional, which shows in particular that W{xi) = Ai > 



(3.7) 



Ai/(xi) + ... + A^/(xw) 



for all Xi ^V. 



□ 
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The coefBcient matrices of the three-term relation can be used to define 

the analog of Jacobi matrices, 



Al. Bi. 



O 



o 



i Bn-2.i ^n-2,i 
AT R 
^n-2,i ^n-l,i 



Ki <d. 



According to the three-term relation, Ji is the matrix representation of the operator 

Pfc mod /(V^), 1 < A: < n, it follows that these 
matrices commute. In other words, we have the following proposition. 

Proposition 3.9. If n is finite, then the Jacobi matrices commute, JiJj = JjJi, 
1 < i.] < d. 

If n is infinity, then the matrices become infinite and we need to consider them 
as operators. Still, the Jacobi matrices formally commute. In the case that A = Nq, 
see 112 . 



4. Examples 

Some examples have been given in the previous sections. To make them more 
concrete, let us mention two classical discrete polynomials. See, for example, jSIISI- 

The Hahn polynomials, Q{x; a, b, N), are discrete orthogonal polynomials defined 
on the set ^ = {0, 1, . . . , N} and are orthogonal with respect to the hypergeometric 
distribution (a + l)j;{b + 1)n-x/ {xl{N - x)!). 



N 

E 



X + a\ f N ~ y + b 



Qnix; a, b, N)Q.m{x; a, b, N) 

_ (-l)"n!(b + l)„(n + a + b+ 1)n+i ^ 
iV!(2n + a + b+ l)(-7V)„(a + 1)„ ' 
Their explicit formulas are given in terms of 3F2 series, 

-n, n + a + b + I, —x 
a+l,-N 



Qn{x;a,b,N) 



3F2 



1 



n = 0, 1, 



n, m < N. 



.,N, 



where 3F2 is defined as the usual 3F2 with the summation terminating at N. The 
Hahn polynomials, Qnix) = Qn{x]a,b, N), satisfy the three-term relation 



-xQn{x) = A„Qn+l{x) - {A„ + Cn)Qn{x) + CnQn-lix), 



where 



A„ = 



{n + a + b+ l)(n + a+ 1){N - n) 



n{n + b){n + a + b + N + l) 



{2n + a + b){2n + a + b+1) ' (2n + a + b){2n + a + b + I)' 

The Meixner polynomials, Mn{x;b,c), are discrete orthogonal polynomials de- 
fined on the set F = No and are orthogonal with respect to the negative binomial 
distribution {b)xC^ /x\, 
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where (a)m denote the Pochhammer symbol {a)m = a{a + 1) ... (a + to — 1). Their 
expUcit formula is given in terms of 2F1 series, 

M„(x) -.^ Mn{x;b,c) = 2-Fi(-n, -x; 6; 1 - c-^), ?i = 0,1,2,... 

which satisfies the three-term relation 

(c - l)xMn{x) = c{n + b)Mn+i - {n + {n + b)c)M,i{x) + nMn-i{x). 

Example 4.1 As a special case of Example 3.1, we have the product Hahn polyno- 
mials of two variables, Qn{x;ai,bi, N)Q,n{x;a2,b2, M), 0<n<N,0<m< M, 
which are orthogonal on the set V — {0, 1, . . . , N} x {0, 1, . . . , M} with respect to 
the weight function 

^x + ai\ /N ~ X + ai\ /y + b\ /N - ij + b^ 
X )[ N-x )\ y )\ N^y 

and the product of Hahn and Meixner polynomials, Qn{x;ai,a2, N)Mm{x;b,c), 
< n < N,m > 0, which are orthogonal on the set V = {0, 1, . . . , N} x No with 
respect to the weight function 

X + ai\ f N ~ x + ai\ {b)y y 



W{x,y) 



X A N^x 

We also have product Meixner polynomials Mn{x; bi, ci)Mm{x; 62, C2), n > 0, to > 
0, which are orthogonal on the set 1/ = Nq with respect to the weight function 
{bi)x{b2)yC^^^ /{xlyl). In the last case, we have = fc + 1 for all k > 0, as for the 
usual continuous orthogonal polynomials. 

For the product orthogonal polynomials, the coefficient matrices Ak^i in the 
three-term relation are given in Example 3.1*. Note that it is easy to get orthonor- 
mal Hahn and Meixner polynomials (multiply by square root of the normalization 
constant), and the product of the orthonormal polynomials gives orthonormal poly- 
nomials in two variables. □ 

Example 4.2. Let V be the "triangle" point set in Example 2.2. Then the orthog- 
onal polynomials P^, < fc < to exist with = fc -I- 1. This is the case that works 
exactly as in the case of the continuous orthogonal polynomials. As an example, 
let us mention the Hahn polynomials of two variables as defined in Jl^. They are 
given by 

, . M\ f 1^n-^m (O-I + l)n(0-2 + I)m(-A^ + X)„ 

(0-3 -I- l)m(0-2 +(T3+2n+ l)n{-N)m 

X Q„(x;a-i,o-2 + 0-3 + 2n + 1, iV - m)Qm{y;cr2,cT3,N - x), 

and are orthogonal on the set V — {{x, y) G Nq : < x + y < N} with respect to 
the weight function 



W(x,y) = 



X + crA /y + crA f N - X — y + aj, 
^ )\ y J\ N -x-y 



In this case rj, = k+1 and the matrix Ak,i and Ak,2 are of the size (fc+ 1) x (fc + 2), 
just as in the usual continuous case. □ 

For extensions of classical discrete orthogonal polynomials to several variables, 
we refer to 01 |H1 El • One can also extend discrete g-orthogonal polynomials to 
several variables. 
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In the following we consider an example in which V contains 8 points and A is 
given as in the right figure of Figure 1. 

Example 4.3. We set 

V ^ {(-1, -1), (0, -1), (1, -1), (-1, 0), (0, 0), (1, 0), (-1, 1), (-1, 2)}. 

Then R[V] = span{l, x, y, , xy, y^, x^y, y^} as shown in Example 2.2. We use the 
method described in Theorem 3.3 to construct orthogonal polynomials on V with 
respect to the linear functional 

^(/) = I E /(^) 

That is, we compute the matrix Ma = (a;^,a;^), where (/, .g) C{fg), and factor 
it as SDS^. The orthogonal polynomials are given as follows: 

P^ix,y) = l, 

Poix,y) = l + 4x, 

Plix,y)=3 + 12x + 22y, 

P^{x, y) = -26 + X + 35x^ - Ay, 

Pl{x, y) = 3 + 3x + x'^ + 6y + 8xy, 

Pi{x,y) = -20 + 31x-a;^ + lly + eOxy + 51?;^ 

P^{x, y) = -20 + 3x + 27x^ - 45y + Axy + h&x^y - 5y^ 

Pl{x, y) = -Qx + Qx^ ~ 50y - 12a;y + I2x^y ~ iOy^ + 20^/^ 

where P^ is a polynomial of degree k. By construction, these polynomials are mu- 
tually orthogonal and become orthonormal upon multiplying by proper constants. 
The corresponding dimensions of Pq, Pi, P2, F3 are 1, 2, 3, 2. □ 

In the case that y is a set of lattice point and V satisfies H2.1|l . we can take 
K = V . Example 4.1 and 4.2 are examples of such a case. If we take ^ = A in 
Example 4.3, we get orthogonal polynomials (x, y) — Pf{x + 1, y + 1), where P^ 
are those in the example. 

Acknowledgments. The author would like to thank one referee for his careful review 
and helpful suggestions. 
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